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A search is performed for the lepton number violating decay B + — > h~n + fi + , where h~ represents 
a K~ or a n~ , using data from the LHCb detector corresponding to an integrated luminosity of 
36 pb -1 . The decay is forbidden in the Standard Model but allowed in models with a Majorana 
neutrino. No signal is observed in either channel and limits of B{B + — > K~ /i + fi + ) < 5.4 x 10 and 
B(B + — > 7r _ /i + /i + ) < 5.8 x 1CF 8 are set at the 95% confidence level. These improve the previous 
best limits by factors of 40 and 30, respectively. 

PACS numbers: 11.30.Fs, 13. 20. He, 13.35.Hb 



Lepton number is conserved in the Standard Model but those with Majorana neutrinos [1] or left-right symmetric 
can be violated in a range of new physics models such as models with a doubly charged Higgs boson [2] . In this 
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letter a search for lepton number violating decays of the 
type B + — > h~ fi + fi + , where h~ represents a K~ or a 7r~, 
is presented. The inclusion of charge conjugated modes 
is implied throughout. A search for any lepton number 
violating process that mediates the B + — > h~ fi + p + de- 
cay is made. A specific search for B + — > h~ decays 
mediated by an on-shell Majorana neutrino (Fig. 1) is 
also performed. Such decays would give rise to a nar- 
row peak in the invariant mass spectrum of the hadron 
and one of the muons [3], m v — m^, if the mass of the 
neutrino is between mjfjrj+ra^ and to_b~ "V- The previ- 
ous best experimental limit on the B + — > K~ (ir~)n + /j, + 
branching fraction is 1.8(1.2) x 10~ 6 at 90% confidence 
level (CL) [4]. 

The search for B + — > h~p + p + is carried out with data 
from the LHCb experiment at the Large Hadron Collider 
at CERN. The data corresponds to 36 pb -1 of integrated 
luminosity of proton-proton collisions at y/s = 7TeV 
collected in 2010. The LHCb detector is a single-arm 
spectrometer designed to study 6-hadron decays with an 
acceptance for charged tracks with pseudorapidity be- 
tween 2 and 5. Primary proton-proton vertices (PVs), 
and secondary B vertices are identified in a silicon strip 
vertex detector. Tracks from charged particles are re- 
constructed by the vertex detector and a set of tracking 
stations. The curvature of the tracks in a dipole magnetic 
field allows momenta to be determined with a precision 
of Sp/p — 0.35-0.5%. Two Ring Imaging CHerenkov 
(RICH) detectors allow kaons to be separated from pi- 
ons/muons over a momentum range 2 < p < lOOGeV/c. 
Muons with momentum above 3 GeV/c are identified on 
the basis of the number of hits left in detectors inter- 
leaved with an iron muon filter. Further details about 
the LHCb detector can be found in Ref. [5]. 

The search for B + — > h~fi + /j, + decays is based on 
the selection of B + — » h ± p, + fi T candidates. The B + — » 
J/ipK + decay with J/0 — > is included in the same 

selection. It is subsequently used as a normalisation 
mode when setting a limit on the branching fraction of 
the B + — > h~ fi + p + decays. The selection is designed to 
minimise and control the difference between decays with 
same- and opposite-sign muons and thus cancel most of 




FIG. 1. s-channel diagram for B — > K~fi (J, 
(B + — > n~fj, + n + ) where the decay is mediated by an on-shell 
Majorana neutrino. 



the systematic uncertainty from the normalisation. The 
only differences in efficiency between the signal and nor- 
malisation channels are due to the decay kinematics and 
the presence of a same-sign muon pair, rather than an 
opposite-sign pair, in the final state. 

In the trigger, the B + — s- /i ± /i + /i =F candidates are re- 
quired to pass the initial hardware trigger based on the 
of one of the muons. In the subsequent software trig- 
ger, one of the muons is required to have a large impact 
parameter (IP) with respect to all the PVs in the event 
and to pass requirements on the quality of the track fit 
and the compatibility of the candidate with the muon 
hypothesis. Finally, the muon candidate combined with 
another track is required to form a vertex displaced from 
the PVs. 

Further event selection is applied offline on fully recon- 
structed B decay candidates. The selection is designed 
to reduce combinatorial backgrounds, where not all the 
selected tracks come from the same decay vertex; and 
peaking backgrounds, where a single decay is selected 
but with some of the particle types misidentificd. The 
combinatorial background is smoothly distributed in the 
reconstructed £>-candidate mass and the level of back- 
ground is assessed from the sidebands around the signal 
window. Peaking backgrounds from B decays to hadronic 
final states, final states with a J/0 and semileptonic final 
states are also considered. 

Proxies are used in the optimisation of the selection for 
both the signal and the background to avoid a selection 
bias. The B + — > J/tpK + decay is used as a proxy for the 
signal. The background proxy comprises opposite-sign 
B + — > h + n + pT candidates with an invariant mass in the 
upper mass sideband and with muon pairs incompatible 
with a J/0 or a ip(2S) hypothesis. The bias introduced 
by using B + — > J/0 if + for both optimisation and as 
a normalisation mode is insignificant due to the large 
number of candidates. 

The combinatorial background is reduced by requiring 
that the decay products of the B have pt > 800 MeV/c. 
Tracks are selected which are incompatible with origi- 
nating from a PV based on the x 2 of the tracks' impact 
parameters (x 2 P > 45). The direction of the candidate 
B + momentum is required to be within 8mrad of the 
reconstructed B + line of flight. The B + vertex is also 
required to be of good quality (x 2 < 12 for three degrees 
of freedom) and significantly displaced from the PV (x 2 
of vertex separation larger than 144). 

The selection uses a range of particle identifica- 
tion (PID) criteria, based on information from the RICH 
and muon detectors, to ensure the hadron and the muons 
are correctly identified. For example, DLL^ T is the dif- 
ference in log-likelihoods between the K and tt hypothe- 
ses. For the B + — > K~ p + p + final state, DLL^ > 1 is re- 
quired to select kaon candidates. For the B + — >• 7r~ /i + p + 
final state the selection criterion is mirrored to select 
pions with DLLrv < -1. The B + -> K~n + p + and 
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B + — » 7r~/i + /i + selections are otherwise identical. 

To reject background events where two tracks that are 
close together in the tracking system share hits in the 
muon detector, a requirement is made on the maximum 
number of muon system hits that two candidate muons 
may have in common. This requirement can introduce a 
bias in the relative efficiency between signal and normali- 
sation channels, as both tracks from the same-sign muon 
pair for B + — > h~ [V y« + will curve in the same direction 
in the dipole field. Simulated events give an estimate of 
0.3% for the effect on the relative efficiency between the 
signal and normalisation channel. In order to avoid se- 
lecting a muon as the pion or kaon, the candidate hadron 
is also required to be within the acceptance of the muon 
system but not have a track segment there. 

After the application of the above criteria the combi- 
natorial background is completely dominated by candi- 
dates with two real muons, rather than by hadrons mis- 
identified as muons. 

The invariant mass distribution and the relevant mis- 
identification rates are required in order to evaluate the 
peaking background. These are evaluated, respectively, 
from a full simulation using Pythia [6] followed by 
Geant4 [7], and from control channels which provide 
an unambiguous and pure source of particles of known 
type. The control channel events are selected to have the 
same kinematics as the signal decay, without the appli- 
cation of any PID criteria. D* — > D°7r, D° — >• Kit decays 
give pure sources of pions and kaons. A pure source of 
muons is selected by using a tag-and-probe approach with 
J/ip—> decays [8]. 

Under the B + — > K~[i + fi + hypothesis, any crossfeed 
from B + — > J/ij)K + decays would peak strongly in the 
signal mass region. The K — > \i mis-id rate is evaluated 
from the above D* sample and the /i — > K mis-id rate 
from the J/ip sample. The later mis-id rate is consis- 
tent with zero. The number of B + —> J/tpK + events ex- 
pected in the signal region is therefore also zero but with 
a large uncertainty which dominates the error on the to- 
tal exclusive background expected in the signal region. 
The B + — > ir~ir + K + decay contributes the most to the 
peaking background with an expected (1.7 ± 0.1) x 10~ 3 
candidates, followed by the B + — » K~ir + K + decay with 
(6.1 ± 0.8) x 10~ 4 candidates. The total peaking back- 
ground expected in the B + — > K~ signal region is 
(3.4lJ 4 2 °) x 10 -3 events with the asymmetric error caused 
by the zero expectation from the B + — > J/ipK + decay. 

Under the B + — > 7r~/.t + /j + hypothesis, B + — > J/ipK + 
decays are reconstructed with invariant masses below 
the nominal B + mass, in the lower mass sideband. 
The dominant background decay in this case is B + — » 
7r~7r + 7r + , where the two same-sign pions are misidenti- 



fied as muons. The B + — > 
level is (2.9 ± 0.6) x 10~ 2 
In Fig. 2(a), the 
distribution for B + 



7T 

events 



peaking background 



l m p+^- ~ to j/v> I < 50MeV/c 2 is shown, after the 
application of the selection. In the B + — > J/ipK + 
sample, there are no events containing more than one 
candidate. An unbinned maximum likelihood fit to the 
B + — > J/ip K + mass peak is made with a Crystal Ball [9] 
function which accounts for the radiative tail. The 
combinatorial background is assumed to be flat, and the 
partially reconstructed events in the lower mass sideband 
are fitted with a Gaussian distribution. The signal peak 
has a Gaussian component of width 20MeV/c 2 , and a 
signal mass window of 5280 ± 40 MeV/c 2 is chosen. The 
B + — > J/ijjK + peak contains 3407 ± 59 signal events 
within the signal window. B + — > J/tjj^ candidates 
were also examined and, accounting for a shoulder in 
the mass distribution from B + — > J/*ipK + , the yield 
observed agrees with the expectation. 

The invariant mass distribution for events 

with — mj/^ ^(2S)| > 70MeV/c 2 is shown in 

Fig. 2(b). Using the same fit model, with all shape pa- 
rameters fixed to those from the above fit, the signal 
peak was determined to contain 27 ± 5 events from the 
B + — > K + n + ii~ decay. The ratio of branching fractions 
between B + ->• J/ipK + and B+ -> K+fi + ^ decays [10] 
and the trigger efficiency ratio predicted by the simu- 
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FIG. 2. Invariant mass distribution of K + n + events after 
the application of the selection criteria. In (a) requiring the 
muon pair to be compatible with coming from a J/V> decay 
and in (b) excluding invariant mass windows around the J/tj) 
and ip(2S) for the muon pair. The curve is the fit to data as 
described in the text. 
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lation, give an expectation of 29 ± 4 B + — > K + fi + [i~ 
decays. The observed yield is consistent with the expec- 
tation showing that the selection does not favour candi- 
dates with a dimuon mass close to the J/tp mass. 

The difference in efficiency between the signal and nor- 
malisation channels was evaluated using Monte Carlo 
simulation samples. The relative selection efficiency 
across the phase space is shown for B + — > K~ in 
Fig. 3. The efficiency of the signal selection in a given 
phase space bin is divided by the average efficiency of 
B + — > J/i/jK + , to yield the relative efficiency for that 
bin. The D* control channel is used to determine the 
PID efficiencies required to normalise B + — »■ 7r~/i+^ + to 
B+^J/iPK+. 

Assuming a signal that is uniformly distributed in 
phase space, the relative efficiency of B + — > K~ \i + n + 
and B + -> J/ipK + was calculated to be 89.1±0.4 (stat)± 
0.3 (syst)%. The relative efficiency of B + — > n~ and 
B+ -> J/i})K + was calculated to be 82.7 ± 0.6 (stat) ± 
0.8 (syst)%. The systematic uncertainties associated 
with these estimates are detailed below. These relative 
efficiencies together with the number of events observed 
in the normalisation channel, give single event sensitiv- 
ities of 2.0 x 10~ 8 (2.1 x 10~ 8 ) in the B+ -> 
(B + — > 7r _ /i + /i + ) case. 

In order to compute the efficiency under a given Majo- 
rana neutrino mass hypothesis, a model for the variation 
of efficiency with m/^ is required. For a given value of 
rrihfj, this is obtained by varying the polarisation of the 
Majorana neutrino in the decay and taking the lowest 
(most conservative) value of the efficiency. 

The dominant systematic uncertainty (under the as- 
sumption of a flat phase-space distribution) for the single 
event sensitivity is the 3.4% uncertainty on the B + — > 
J/tpK + branching fraction. The statistical uncertainty 
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FIG. 3. Relative efficiency between the B + — ¥ K~ /i + /i + signal 
and the B + — > J/i[>K + normalisation channel. The plot has 
been symmetrised over the diagonal. 



TABLE I. Sources of systematic error and their fractional 
uncertainty on the relative efficiency. 



Source 


B+^ K~ij, + h + 


B + ^tx-^+ 1 j+ 


B(B+^ J/4>K+) 


3.4% 


3.4% 


J/tpK + yield 


1.7% 


1.7% 


B+-+ J/tpK + fit models 


1.6% 


1.6% 


Simulation statistics 


0.4% 


0.6% 


IP modelling 


0.2% 


0.2% 


PID modelling 


0.1% 


0.8% 


Trigger efficiency 


0.1% 


0.1% 


Tracking efficiency 


0.1% 


0.1% 



on the B + — > J/ij)K + yield gives an additional system- 
atic uncertainty of 1.7% and the uncertainty from the 
model used to fit the data is 1.6%. The latter is eval- 
uated by changing the Crystal Ball signal function used 
in the fit to a Gaussian and the polynomial background 
function to an exponential. 

There are several sources of uncertainty associated 
with the calculation of the relative efficiency between the 
signal and normalisation channels. In addition to the 
statistical uncertainty of the simulation samples, there 
are systematic uncertainties from: the differences in the 
effect of the IP selection criteria between the simula- 
tion and data; the statistical uncertainty on the mea- 
sured PID efficiencies; the uncertainties associated with 
the simulation of the trigger; and the uncertainty in the 
tracking efficiency. In each case the systematic uncer- 
tainty is estimated by varying the relevant criteria at the 
level of the expected effect and re-evaluating the relative 
efficiency. For the B + — » ir~fi + fi + decay, there is an ad- 
ditional uncertainty from the correction for the relative 
kaon- and pion-identification efficiencies. The systematic 
uncertainties averaged over the three-body phase space 
are given in Table I. 

A limit on the branching fraction of each of the B + — > 
h~ fi + fjb + decays is set by counting the number of ob- 
served events in the mass windows, and using the sin- 
gle event sensitivity. The probability is modelled with 
a Poisson distribution where the mean has contributions 
from a potential signal, the combinatorial and peaking 
backgrounds. The combinatorial background is uncon- 
strained by measurements from the simulation or the 
opposite-sign data. The number of events in the upper 
mass sideband is therefore used to constrain the con- 
tribution of the combinatorial background to the Pois- 
son mean. The upper mass sideband is restricted to 
masses above m/j^ > 5.4 GeV/c 2 such that any peak- 
ing background component can be ignored. In both the 
B + — > K~fi + fi + and B + — > 7r~/i + ^t + cases no events are 
found in either the upper or lower mass sidebands. This 
is consistent with the observation of three opposite-sign 
candidates seen in the B + — > K + (i + [i~ upper mass side- 
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band (Fig. 2) and two candidates in the B + — » 7r + /x + /x~ 
upper mass sideband. The peaking background estimates 
are explicitly split into two components, the contribution 
from B + — > h~h + h + decays and that from B + — > J/ipK + 
decays. The latter has a large uncertainty. The cen- 
tral values for both peaking background components are 
taken from the estimates described above. 

Systematic uncertainties on the peaking background, 
single event sensitivity and signal-to-sideband scale fac- 
tor are included in the limit-setting procedure using a 
Bayesian approach. The unknown parameter is inte- 
grated over and included in the probability to observe 
a given number of events in the signal and upper mass 
window. 

In the signal mass windows of B + — > K~fi + fi + and 
B + — > 7r~/i + /i + no events are observed. This corresponds 
to limits on the B + —> branching fractions of 

B{B+ -> K~ /!+ fi+) < 5.4(4.1) x 1CT 8 at 95% (90%) CL, 
B{B+ -> 7T - /i + m + ) < 5.8(4.4) x 10~ 8 at 95% (90%) CL. 



cay and, using B + — > J/ijjK + as a normalisation chan- 
nel, the present best limits on B(B + — » K~/j + /.i + ) and 
B(B + — > 7r~/i + /i + ) are improved by factors of 40 and 30, 
respectively [4]. 
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FIG. 4. The 95 % CL branching fraction limits for B + -> 
K~ fi + n + (light-coloured line) and B + — > -k~ fi + fi + (dark- 
coloured line) as a function of the Majorana neutrino mass 



The observation of no candidates in the sidebands as 
well as the signal region is compatible with a background- 
only hypothesis. The dependence of the limit in 
models where the Majorana neutrino can be produced 
on mass shell is shown in Fig. 4. The shapes of the limits 
arise from the changing efficiency as a function of mass. 

In summary, a search for the B + — > K~ fi + [i + and 
B + — > 7r~/i + /i + decays modes has been performed with 
36 pb -1 of integrated luminosity collected with the LHCb 
detector in 2010. No signal is observed in either de- 
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